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Abstract 



> 

C^ ' We present a fully nonequilibrium calculation of the low temperature trans- 

^2^ ' port properties of a quantum dot in the Kondo regime when an AC potential 

O 

Q ' is applied to the gate voltage. We solve a time dependent Anderson model 



with finite on-site Coulomb interaction. The interaction self-energy is calcu- 
lated up to second order in perturbation theory in the on-site interaction, in 



r^ ' the context of the Keldysh non-equilibrium technique, and the effect of the 

o' 

J-^ ' AC voltage is taken into account exactly for all ranges of AC frequencies and 

Lj AC intensities. The obtained linear conductance and time- averaged density 

C^ ' of states of the quantum dot evolve in a non trivial way as a function of the 

AC frequency and AC intensity of the harmonic modulation.. 
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I. INTRODUCTION 

Recent experiments [|l]-0| showing Kondo behavior in the low temperature transport of 
quantum dots (QD's) have opened a new arena for the study of strongly correlated electrons 
in artificial systems. Kondo effect in dilute magnetic alloys appears as a crossover from 
weak to strong coupling between delocalized electrons of the host non-magnetic metal and 
the unpaired localized electron of the magnetic impurity as the temperature (T) is reduced 
well below the Kondo temperature (Tx) @]. This crossover leads to the formation of a 
singlet state between the unpaired localized electron in the impurity and electrons in the 
host metal. 

It is important, however, to emphasize the main differences of the Kondo physics in QD's 
with respect to bulk magnetic impurities. The parameters which define the T^- in QD's can 
be changed in a controlled way by applying the appropriate combination of gate voltages. So, 
it is possible to study either Kondo or Mixed Valence regimes in the same sample. For this 
to be possible, there is an important requirement: the charging energy and level separation 
of the QD must be significantly larger than the level broadening due to the coupling to the 
leads. More importantly, the study of Kondo physics in QD's opens a new road to the study 
of nonequilihrium many-body phenomena, a relatively young and rich area in contemporary 
condensed matter physics. 

In this paper, we address the issue of a QD driven out of equilibrium by means of an AC 
voltage. More specifically, we study theoretically the low temperature transport properties 
of a QD with an AC voltage applied to the central gate. We use a time dependent version 
of the Anderson model. In its simplest formulation, the Anderson model, valid for both 
Kondo and Mixed Valence regimes in bulk systems, describes a single discrete level with on- 
site electron-electron interaction coupled to a band. The model describes different physical 
regimes which, for QD's, are determined by the following parameters: (i) the energy differ- 
ence between a discrete level in the QD (eo) and the Fermi energy of the leads {Ep)- (ii) the 
tunneling coupling (F) between the discrete level in the QD and the electronic states in the 



reservoirs, (iii) the QD charging energy (on-site interaction U), i.e., the energy necessary to 
add an electron to the QD. The relevant energy scale is Tk ^ V2t/re"''[l^^"'ol(^+'o)l/^'^^, 
which is related to the binding energy of the many body state ||^. For T < Tk-, the Kondo 
regime is reached when eo < Ep — F and eo > Ep — f/ + F and the Mixed Valence regime 
is established for Ep — F < cq < Ep and Ep — f/ < cq < Ep — U + V. In the Kondo 
regime at T = 0, the low energy excitations (quasi-particles) produce a peak at Ep (Kondo 
resonance or Abrikosov-Suhl resonance) in the density of states (DOS) f^. One electron 
at Ep becomes scattered by the QD undergoing a phase shift which is proportional to the 
displaced charge 5n (Friedel-Langreth sum rule (FLSR) 0) and the linear conductance for 
a QD symmetrically coupled to the leads takes the value Q = (2e^//i)sm^(7r5n). In the 
symmetric case {eo = ^), Sn = 0.5 leads to a perfect transparency of the QD. For any 
chemical potential between eo and eo + U the QD has the linear conductance as a function 
of the chemical potential of an almost perfectly open channel 2e^ /h [|]-§]. However, as 
temperature increases, inelastic scattering processes reduce the DOS at Ep (i.e., the linear 
conductance) at the Kondo valley and eventually two peaks at eo and t^ + U appear for 
T » Tk. 

As we already mentioned, new questions arise when driving the QD out of equilibrium 
|9|^T|. When this is done by means of the application of a finite DC voltage bias, the linear 
conductance is reduced and the Kondo peak in the DOS splits ITIHT^. More sophisticated 



configurations of QD's in the Kondo regime constitute a growing area of intense investiga- 
tions, both from the theoretical and experimental sides. Time dependent Kondo physics 
1^ , Kondo physics in integer-spin QD's |jl5|, QD's embedded in Aharonov-Bohm rings |1^ 



or double QD systems [0 , are examples of such configurations. 

We focus on the study of the transport properties of an Anderson Hamiltonian with a 
time dependent resonant level eo(t) = eQ + VAccos{ujQt). This can be achieved experimentally 
by means of a time dependent central gate voltage capacitively coupled to the QD. In the 
high temperature regime these type of experiments have indeed been carried out leading 



to the observation of photon-assisted processes in the Coulomb blockade regime ||18|,|T9 



This regime has been studied theoretically as well pO|. In the same way, there has been 
some theoretical effort devoted to studying the AC transport at very low (T « T^) and 
intermediate temperatures (T ^ Tk) in the Kondo and Mixed Valence regimes pl|-p7|. 
Recently transport in an AC driven quantum dot at low temperatures has been measured 
as well |2§. 

In this work we clarify the role of an AC voltage in the Kondo effect in QD's. We 
concentrate on the low temperature regime so that a Fermi liquid theory is adequate in 
spite of the fact that it overestimates the width of the Kondo peak, i.e., the Tk ^^^- In 
fact our Fermi liquid approach gives a resonance width which decreases algebraically with U, 
instead of having an exponential decay as given by scaling calculations [Q . The occupation of 
the QD as well as all the relevant quantities in transport have to be calculated by using non- 
equilibrium propagators. Finally, some approximation is needed for calculating the Green's 
function of the QD. We use a finite-U perturbation theory for the Green's function of the 
impurity. The Anderson Hamiltonian is solvable through the Bethe Ansatz [^, Numerical 
Renormalization Group |]^ and Quantum Monte Carlo methods [^ but a reliable and 



simple method to obtain dynamical properties at low temperatures in the whole range of 
interactions (U/T) is not available. Previous efforts have concentrated on the U ^ oo 
limit where a Non-Crossing Approximation (NCA) |TI|Jl^ , |2l|j2^ can be made for high and 



intermediate temperatures. However, such approximations do not give the exact local Fermi- 
liquid properties as T — i> 0. They cannot describe the transition from the weak-correlation 
to the strong-correlation regime. On the other hand, finite-f/ perturbation theory |9|, |33|j3^ 
describes the symmetric case properly, but presents clear anomalies away from this special 



situation (which can be overcome by means of an interpolative self-energy [p!0| , p5| ). In this 



paper we will restrict ourselves to the symmetric case which is the relevant one for the 



experimental information available ||28[ . The study of the asymmetric case will be analyzed 
elsewhere. 

The main difficulty for our purpose resides in the determination of the QD Green's 



function, and specifically of its self-energy. In a previous paper p5| we proposed an ansatz 



for the modification of the QD Green's function due to an harmonic modulation. Here, we 
improve our previous description, vahd in the hmit of small interaction U, and extend the 
calculation to finite temperatures. 

The paper is organized as follows: in Sec. II we describe the theoretical model and de- 
duce the expressions for the self-energy and the time-averaged spectral density. In Sec. Ill 
we present the results for the time-averaged density of states and the linear conductance at 
finite temperatures and for different AC frequencies and AC intensities of the harmonic mod- 
ulation. Moreover, we compare with previous experiments and theoretical models. Finally, 
Sec. IV summarizes the main conclusions of the paper. 

II. THEORETICAL MODEL 

A. Keldysh Formalism 

The application of a time dependent component to the energy level in the Anderson 
Hamiltonian (See Eq. 6, below) breaks the time translational invariance of the system 
and, then, we need an approach capable of addressing this fully non-equilibrium situation. 
When the time dependent perturbation acts for a while, the system does not recover its 
thermodynamic equilibrium after the perturbation is over. The whole process has not the 
symmetry between t ^ — oo and t —>■ oo and, then, it is not possible an equilibrium expansion 
in terms of expectation values. Nevertheless, the problem can be solved by allowing the 
system to evolve from — oo to the moment of interest i and then continuously evolve from 
t = i back to t — i> — oo. In this way all the expectation values are evaluated in a well- 
defined state which was prepared in a remote past. This special complex time-contour is 
the main ingredient of the non-equilibrium Keldysh formalism [^ . By using this complex 
time-contour (see Fig. |l|) the formal expression for the Green's function relative to a total 
Hamiltonian H = Hq + Vint which includes an interaction potential Vint is given by 

G,,.(t,t') = (T,[5_(-oo,oo)5+(oo,-oo)d.(t),4(tO]), (1) 



where S'+(oo,— oo) is the usual S-matrix, exp[—i J^^Vint{T)dT] and S''''(oo, — oo) = 
S'(— 00,00). d^{t), (i^(t') are operators in the interaction picture, ^^(oo,— 00) depicts the 
causal branch (from —00 to 00), while S'_(— 00, 00) depicts the anticausal branch (from 00 
to —00). The Green's function defined in Eq. (^ is a matrix 2x2, since both times, t and 
t' can belong to the causal branch (causal Green's function, G^^{t,t') ) or to the anticausal 
one (anticausal Green's function, G^^{t, t')) and it is possible to have one time in the causal 
branch and the other one in the anticausal branch (which defines the lesser, G'^^(t,t') and 
greater G^^{t,t') Green's functions). Tc is the contour-ordering operator which arranges 
the operators on the closed time-path in this way: operators with time labels later on the 
contour are moved left of operators of earlier time labels. In this way, the non-equilibrium 
perturbation theory has the same structure as the equilibrium perturbation theory keeping 
in mind the label of the times. The 2x2 matrix for the Green's function reads 



Gda(t,t ) = 

where 



^G^,(t,t') G< (t,t')^ 



G>(t,t') G^(t,t' 



(2) 



Gj,(t,t') = -z{T{d^{t),di{t')}), G%{t.t') = -2(T{rf.(t),4(t')}), 

G<,(t,0 ^ i{dl{t')d„{t)), G>^{t,t') ^ -z{d^{t)dUt')). (3) 

Here, T is the causal time-ordering operator, and T is the anticausal time-ordering operator. 
In the conventional Keldysh matrix formulation of the perturbation theory, one does not 
work directly with the Green's function defined on the complex time-contour, but with a 
linear combination of these four possible time orders. The usual linear combinations are 
(similar relations hold for the self-energies) 

Q^(t,t') = e{t - t')[G>^{t,t') - G<^(t,t')], 

GlAt,t') = e{t' - t)\GlXtX) - Gl„{t,f)l (4) 

where G\^{t^ t') is the retarded Green's function, G'^^it, f) is the advanced Green's function, 
and G^^(t,t'), G^^(t,t') are the so-called lesser and greater Green's functions, respectively. 



B. Hamiltonian 

The time dependent Anderson Hamiltonian is: 

H = Hieads + Hqd + Hsd + HAc{t), (5) 

where 

k£{L,R},cT 

Hqd = "^eo^adlda + Ud\d^dldi, 

(7 

Hsd = J2 ^kci^ada + V;*4cfc,a, 

k&{L,R},a 

HAc{t) = J2 VAccosuot dld„. (6) 

Vac and ojq are the AC intensity and AC frequency of the AC potential respectively. d\ 
creates an electron with spin a in the QD, while c|, ^ creates it in the lead with energy e^ [k 
labels the rest of quantum numbers). The AC voltage modulates in time the relative position 
of the QD level eo,o- with respect to Ep- An eventual breakdown of the spin degeneracy would 
be represented by eo,cr 7^ £0,^. The coupling Vk between the QD and the leads produces a 
broadening r^(^)(e) = -2Im[^^J/^\e + i5)] = 27i EkeL{R) {Vk^^ie - ek), where S^^^ is 
the hybridization single-particle self-energy. Hereafter, for simplicity, we consider the wide 
band (WB) limit approximation which neglects the principal value of the hybridization 
self-energy and considers the imaginary part to be an energy independent constant, i.e., 

C. Model 

Here, we discuss the procedure for obtaining the QD Green's functions which allows us 
to obtain the spectral density of the QD and the linear conductance. 

In the remote past, the QD is decoupled from the leads. The coupling between different 
regions (the contacts and the central region) is treated as a perturbation by means of stan- 
dard equilibrium perturbation theory. In a first step, the effect of the on-site interaction is 

7 



included via a Hartree mean- field approximation (see Appendix A). The time modulation of 
the QD level is treated via non-equilibrium perturbation theory, since the time translational 
invariance is broken by the AC voltage. At this point, we include the correlation effects 
by computing the on-site interaction self-energy (lesser and greater) up to second order by 
means of the diagrams of Fig. 0. These diagrams are evaluated by using the previous lesser 
and greater Green's functions as bare propagators (Appendix B). These bare propagators 
include the coupling between the QD and the contacts, the time dependence of the QD level 
and the on-site interaction in the Hartree approximation. Once the correlation self-energy 
has been calculated, the QD retarded Green's function is obtained by means of the time 
dependent Dyson equation (Eq. 13). Finally, the time-averaged spectral density (Eq. 21) 
and the linear conductance (Eq. 22) are calculated from the QD retarded Green's function. 

D. Correlation self-energy 

The starting point for the derivation of the correlation self-energy in the presence of 
the AC potential, consists of calculating the lesser and greater QD Green's functions, in 
the Hartree approximation, including coupling to the leads. In the absence of the time 
dependent potential, the retarded and advanced QD Green's functions have the following 
expressions (see Appendix A) 

g2;(t - f) = TiOi±t T t') exp-'/*'^*^(^""^+^"^-^'^-",Hr./2) ^ ^^^ 

^d,(T = {dl.it)do-(t)) being the QD occupation. Note that these QD Green's functions 
have been calculated taking into account the coupling self-energy (which is given by 
TiY^a€L,R^a/'^) and the interaction in the Hartree approximation (given by Uud,^). Now, 
if one also considers a time modulation of the QD level, the retarded and advanced QD 
Green's functions have the following forms (see Appendix B) 

GZit, t') = e-^^(--oi-sm.oi')gr,a (^ _ ^,y (g^ 



The lesser and greater QD Green's functions can be obtained through the well known relation 



Glait.t') = J dh J dt2Gl^{t,U)J:f^{tuh)Gl^{h,t'), (9) 

where T,ff'(ti,t2) are the lesser and greater coupling self-energies defined in Appendix A 
(Eqs. A2 and A3, respectively). Now we include correlation effects up to second order in the 
on-site Coulomb interaction (see Fig. |^). The new lesser and greater correlation self-energies 
are calculated by means of the diagrams of Fig. ^ with bare lines which are given by the 
propagators of Eq. @ (analytical expressions are given in Appendix B, Eqs. B11-B14): 

^lf\t,t') = ^U'G>At,t')G<,{t',t)G>At,t^ (10) 

and 

^lfit,t') = -zU'G<Jt,t')G>,{t',t)G<,{t,t'). (11) 

The retarded self-energy 

^'S^ (t, t') = e{t- 1') [s<f ) {t, t') - E>f {t, t% (12) 

is obtained from Eqs. ([10|) and ([iTl). 

E. Dyson equation 

The next step for deriving the retarded QD Green's function consists of solving the 



retarded time dependent Dyson equation. By using Eq. ([T^) for the self-energy, one can 
write 



.d 
I 



- (-^oAt) - ^ E rj2)]G^f (t,tO = 5it - t') + /rftiS^f (t,tOG2,!f)(ti,tO, (13) 



aGL,R 



where eo,o-(t)=eo,o- + Urid^^it) + VAccosujQt. Eq. (|T^) simplifies considerably by making the 
gauge transformation 



9 



In the presence of time modulation, the retarded Dyson equation becomes 

— g„(t, t') = - J^^ dt,K^{t, h)g^{h,t'), (15) 

which is defined only for t > f due to the 9 function appearing in Eq. (|14D. Ka-(t,t') is 
the kernel of the integro-differential time dependent Dyson equation which is related to the 
retarded self-energy through the relation: 

When t = t', an additional condition must be imposed in Eq. (|T3|): G^^\^ {t,t) = 
—i{{dcr(t),d'l{t)}) = —i, where {} is the anticonmmutator. This condition implies that 
the solution of Eq. (|I5|) when t = t' is gait, t) = 1. 

We solve Eq. (|TB|) by discretizing the temporal variables, the partial derivative is replaced 
by the finite difference 

_^- /, ,M ga(m,n)-g^(m-l,n) 
where 6 is the grid spacing in time space. The integral is converted into a sum 

rt "I 

- I dtiKait, ti)ga{tu t') ^ -5 E CkKaim, k)ga{k, n). (18) 

•^*' k=n 

Now, the indexes m, n and k replace the time arguments which appear in the (Eq. ^). The 
coefficients Ck are equal to 1 except when m = fc or n = A; in which Ck = \- In this way, the 
discretized time dependent Dyson equation has the following form, 

m 

ga{m, n) = gaim - l,n) - 5^ Yl CkKa{m, k)ga{k, n). (19) 

k=n 



Eq. (|T9D constitutes a set of linear equations that can be solved by standard numerical 
techniques. Its solution gives the retarded QD Green's function which is used to study the 
transport properties of the system in the next section. 

10 



F. Time-averaged spectral density and Linear Conductance 

The time dependent spectral density Pa{^, i), being i = ^i^, is defined as the imaginary- 
part of the Fourier transform with respect to r = t — t' of the retarded QD Green's function 

p,(e,t) = —Im / G''£\i+ T/2,i- r/2)e*^^rfr. (20) 

Since the measurement of the hnear conductance imphes a time-average in i, we work with 
the time-averaged spectral density which reads 

{A^{e)) = ^f^dipAe,i). (21) 

zvr JO 

The linear conductance [^ at finite temperature, in terms of the time-averaged spectral 
density, is given by 

where /(e) is the Fermi-Dirac distribution function. 



III. RESULTS 

We solve numerically the set of linear equations (0) for a QD in the Kondo regime 
(symmetric case cq = —U/2 with symmetric couplings F/, = F/j = F and U=2.57rF) and a 
temperature T=0.05F ||38| for different AC frequencies and AC intensities. From the solu- 
tion of Eq. (|l^) for the retarded Green's function we obtain the time-averaged density of 
states (Eq. 21) and the linear conductance (Eq. 22). 

The main effect of the AC potential consists in a reduction of the time-averaged DOS at 
Ep- This reduction can be interpreted as decoherence induced by AC excitations, either 
by real photon-assisted induced excitations at large AC frequencies 1^^ or virtual spin-flip 



cotunneling processes at small AC frequencies [^. These processes introduce a quenching 
of the Kondo peak causing a deviation of the linear conductance from the unitary limit. This 
suppression induced by the AC potential depends on two competing mechanisms regulated 

11 



by the AC frequency and by the absorption rate of photons. If the AC frequency is small 
(tuo < Tk) , a negligible reduction will take place even for a large absorption rate of photons 
(the probability of absorption depends on Bessel functions with argument /5 = Vac/^o, as 
shown in appendix B). On the other hand, for large AC frequencies the mechanism of sup- 
pression (AC excitations) becomes ineffective as /? -^ 0. 

To analyze in more detail the behavior of the time-averaged DOS as a function of the AC 
parameters, we plot in Fig. ^ the time-averaged DOS for three different AC frequencies, 
uo = r/4 ^ 2Tk (Fig. |a), uo = r/2 ^ 4Tk (Fig. |b) and ujq = T ^ 8Tk (Fig. |c). For 
each AC frequency we consider different AC intensities corresponding to (3= 0, 0.25, 0.5 and 
0.75. Fig. Ija corresponds to Uq = r/4 ^ 2Tk and different AC intensities Vac ~ Tk/2 
{(3 = 0.25), Vac ^ Tk {(3 = 0.5) and Vac ~ ^Tk/2 {f3 = 0.75). In all these situations the 
Kondo peak is slightly reduced, i. e. the dynamics of the correlated collective state is prac- 
tically not affected by the AC potential since both, the AC frequency and the AC intensities 
are of the order of Tr-. Furthermore, there is no evidence of replicas of the Kondo peak in 
the time-averaged DOS. By doubling uq (Fig. ^), the Kondo peak undergoes a stronger 
reduction as the AC intensity grows. In the present case, the two external energy scales, Uq 
and Vac, are larger than T^ yielding to a stronger modification of the time-averaged DOS 
at Ep. However, a total suppression of the Kondo resonance is only reached for our largest 
AC frequency (Fig. ^, ujq = T). Even for a small value oi (3 = 0.25, the Kondo resonance 
is reduced remarkably. For the highest AC amplitude (/? = 0.75) the Kondo resonance has 
been destroyed completely. Another important effect induced by the AC signal in the time- 
averaged DOS is that the Kondo peak develops satellites. These satellites become apparent 
in Figs. ^D and |^c. In these cases both c^o and Vac are much larger than the relevant 
energy scale of this problem given by T^- In Fig. ^ even the second satellites are resolved 
at /5 = 0.5 [p9|. In a mean field model, one should expect an appreciable increase of the 



satellites of the main peaks in the time-averaged DOS as f3 increases |^,^. However, in the 
present case, the satellites grow very slowly with f3 due to the two competing mechanisms 
previously discussed. It is important to note that, the mean field resonances located at eo 
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and eo + ^ do not display satellites in the time-averaged DOS, since ujq is smaller that their 
widths which are of the order of 2r. 

To illustrate the previous discussion, we plot in Fig. | the conductance as a function of 
Vac for four different AC frequencies. The height of the Kondo peak falls off as the applied 
AC intensity increases. Moreover, the suppression of the time-averaged DOS at Ep occurs 
in a very similar way for the low (cjq = r/4) and intermediate AC frequencies (tuo = r/2, 
uq = T). This decrease is slower for the largest AC frequency (cjq = l.TSF), since for each 
AC amplitude the absorption rate of photons , which depends on /?, is the smallest for the 
largest AC frequency. 

More interesting is the analysis of the conductance as a function of uq, at a fixed AC in- 
tensity (in Fig. ||a for Vac < T and in Fig. |^ for Vac > T). In both cases we find that 
there are two regimes separated by a threshold frequency, uJt, where the conductance has 
a minimum. The presence of this minimum can be understood, again, by the competition 
of two opposite effects. In the adiabatic limit, where cuq — ^ 0, the effect of the applied AC 
potential has a negligible effect on the dynamics of the correlated collective state. In this 
case, the quenching of the Kondo peak is no longer effective and, therefore, the AC only 
induces a small decoherence in the correlated state 1^^. In the very high AC frequency 
limit, where ujq —>■ oo (in this limit (3 —>■ 0), the unpaired electron has a negligible absorption 
probability. Thus, it is obvious from the previous discussion that the conductance versus ujq 
should display a minimal value. For low AC intensities (Fig. ^), the threshold frequency 
depends on the AC amplitude : for Vac = 0.25T Ut ^ O.SF and for Vac = O.SF Ut ~ 0.625r. 
In both cases, for AC frequencies lower than ut the conductance decreases as c^o increases 
whereas this tendency is reversed for AC frequencies larger than ut- Note that for ujq < ojt 
the slowest drop of conductance corresponds to the smallest AC amplitude. For ojq > Ut the 
conductance increases with increasing ujq. Now, the conductance recovers its undriven limit 
for smaller AC frequencies in the case of the smallest AC intensity as one should expect. The 
position of the thresholds frequency depends on the competition of the two effects previously 
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discussed. This competition is illustrated in Fig. Qd, where we plot the time-averaged DOS 
at Ep for three AC frequencies t^o < ujt. By increasing the AC frequency from ujq = 0.375r 
up to ujq = O.Sr the Kondo resonance is suppressed more effectively than by raising the 
AC frequency from c^o = O.SF up to c^o = 0.625r where the reduction of the Kondo peak 
is almost negligible. Fig. ^a depicts the conductance vs. uq for Vac > T. Remarkably, if 
the threshold frequency is normalized to the AC amplitude we find that Ut decreases ex- 
tremely slowly for the three AC intensities {ujq/ Vac — 1). Here, again the suppression of 
the conductance is slower for the lowest AC amplitude. Furthermore, in the region where 
the conductance is an increasing function of the AC frequency, this increase is slower for 
the highest AC intensity. The previous behavior is illustrated in Fig. ^3 where we plot the 
time-averaged DOS for several AC frequencies at the strongest AC intensity Vac = 2r. 



Let us now compare our results with the ones obtained with different schemes p7| , p^ : Kamin- 
ski et al. |^ predict no suppression of the Kondo peak for the symmetric case (eo = —U/2) 
in clear disagreement with both our theory and experimental results [^. In asymmetric 
cases, they predict a monotonous decrease of the height of the Kondo peak as a function of 
uoq (in the low AC frequency regime at fixed Vac)- Their results (for cq 7^ —U/2) are similar 
to our results (for eo = —U/2) for AC frequencies below ojq ~ V /2 and for low AC intensities 
(solid line in Fig. ^) where the Kondo peak is not strongly suppressed. The rest of the 
cases discussed here are away from the range of validity of Ref. p] 



Using a NCA, Nordlander et al. pj] find a nonmonotonous behavior of the conductance vs. 
ujq and the existence of a minimum, in qualitative agreement with us. Caution is needed, 
however, in comparing both works due to the different approximations involved and the 
different regimes of validity for both calculations. Whereas our perturbative calculation is 
valid for the symmetric case, i. e. particle-hole symmetry, and low temperatures, the NCA 
considers the limit U ^>- 00 {i. e. strongly asymmetric case) and temperatures of the order 
and higher than Tk- In that sense, both calculations should be regarded as complementary. 
In Fig. 1^ we plot the time-averaged DOS at Ep as a function of T in the absence of the 
AC potential (solid line) and in the presence of the AC potential for ujq = F, Vac = 0.25F 
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(dashed line) and Vac = O.SF (dashed hne). In the absence of the AC potential, we get the 
well known strong reduction of the DOS at Ep when the temperature increases, as expected 
^. For the lowest AC intensity case (dotted line) a strong temperature dependence is ob- 
served similarly to the case without AC potential (solid line). However, such a dependence 
changes for the strongest AC intensity (dashed line). In this case, the time-averaged DOS 
at Ep has been already so strongly reduced at very low temperatures (T = 0.02r) that an 
increase of temperature produces only a small further reduction of the linear conductance. 
This result can be expected for a system where the Kondo effect is already very weak at low 
temperatures . This behavior of the Kondo peak height as a function of temperature has 
been confirmed experimentally (Fig. 8 in Ref 



IV. SUMMARY 

In conclusion, new features are found in the finite temperature transport (T << T^) 
properties of QD's in the Kondo regime as an oscillatory gate voltage is applied. By solving 
the time dependent Dyson equation we obtain the QD retarded Green's function and the 
time-averaged density of states within the framework of the Fermi liquid theory. The inter- 
action self-energy in our model is calculated, in the context of the Keldysh non-equilibrium 
technique, by perturbation theory up to second order in the on-site interaction and the ef- 
fect of the AC potential is taken into account exactly for all ranges of AC frequencies and 
AC intensities of the AC potential. The Kondo resonance is modified by the external AC 
voltage in a different way depending on the range of AC frequencies studied. We find two 
different AC frequency ranges of opposite behavior, separated by a threshold frequency ujt 
where the linear conductance is minimum, ut depends on the AC intensity and moves to 
higher values as Vac increases. At small AC frequencies, and fixed AC intensity, the Kondo 
peak decreases as c^o grows. Once the threshold frequency is reached the opposite behavior 
is found and the Kondo peak increases as /3 = ^^ -^ 0. Our results quahtatively agree 
with previous theories and complete the range of parameters (arbitrary AC frequencies and 
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AC intensities at finite on-site interaction) and the regime of temperatures (temperatures 
below Tk). We also analyze the behavior of the Kondo peak as a function of temperature at 
fixed AC frequency. For large AC intensities, we obtain a very small decrease of the Kondo 
peak as temperature increases, while for small AC intensities, a much larger quenching of 
the peak is obtained. 
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APPENDIX A: QD GREEN'S FUNCTIONS IN THE ABSENCE OF AC 

POTENTIAL 

In this appendix, the QD Green's functions in the absence of time modulation are ob- 
tained. First of all, the QD Green's function coupled to the leads is calculated, and after- 
wards the interaction in the Hartree approximation is included. A very simple calculation 
yields the exact lesser and greater Green's functions for Hieads, 

g<'^'\t - t') ^ z(4(t')c,(t)) = z/(e,)e-'=(*-*'), 

^>'(°)(t - to ^ -z(c,(t)cl(t')) = -^(1 - /(e,))e-^(*-*'), (Al) 

where /(e^) is the Fermi-Dirac distribution function. The lesser and greater hopping self- 
energies including hopping between the QD and the leads are written in terms of the previous 
Green's functions as [ p^ , 

K,{ti-t2)= E V:g<{h-t,)V, , = zJ: f ^e-'^'^-'^^ X h/n{e)T'^/%) , (A2) 

keL,R L,R'' ^^ 
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^Uti-t2)= E V:g>itl-h)V,,= -^J:|^e-^'^-'^^x{l-h/n{e))T^^^{e). (A3) 



de 
keL,R L,R- 

where r^(e) = 2'Kj2keL |VfcP<^(e — Cfc) and a similar expression is obtained for r'^(e). The 
retarded and advanced self-energies fulfill the relations 

s;,(ti - 12) = 9ih - h)i^Uti - h) - s<,(ti - h)), 

S:,(ti - h) = eih - ti)(S<,(ti - h) - S> (ti - t,)). (A4) 

The Dyson equation for the retarded and advanced QD Green's functions is 

9Z(i - t') = iC^'\t -t') + J dti J dhgX/'\t - U)T.Z{h - h)gZih - t') . (A5) 

Here f?^'" (t — t') are the QD Green's functions for Hqd (Eq. (6)) without the on-site 
repulsion term, 

i,:f\t - t') = T^9i±t T t')e---(*-*'). (A6) 

In the absence of time dependence, it is advantageous to Fourier transform, getting 

C(^) = [(C/°V^-Sr(e)]-\ (A7) 

where 

K',%e) = E J^^'' . = E A.(e) T |r„(e) . (A8) 

k&L,R ^ ^k^^V a£L,R ^ 

In the WB limit approximation the previous self-energies become Ss^(e) = =FZ]aeL,fi fTo,. 
By inserting these self-energies in Eq. ( [A7|) , the retarded and advanced QD Green's functions 
coupled to the leads are 

C(^) = -t\t^ f^- (A9) 

One can Fourier transform back to get the time dependence of the retarded and advanced 
QD Green's functions 
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Now the on-site interaction self-energy in the Hartree approximation is calculated via per- 
turbation theory up to first order in U: 

S2f ^ = Un,,^^ (All) 

where Ud^^ = {da{t)da{t)). By using the Dyson equation it is straightforward to get 

g2^(t,t') = Ti0i±t T t')e~'^''''*^(^'^"^+^"^''^^'^«6^.«^'^/'^ . (A12) 

APPENDIX B: QD GREEN'S FUNCTIONS IN THE PRESENCE OF AC 

POTENTIAL 

In this appendix the analytical expressions for the lesser and greater QD Green's func- 
tions in the presence of an AC potential, coupled to the leads in the Hartree approximation 
are derived. As we said in Sec. II, one needs to obtain these propagators in order to have 
an expression of the interaction self-energy which is given by Eq. ([T2|). However, the lesser 
and greater QD Green's function are given as a function of the retarded and advanced self- 
energies (see Eq. (^). Therefore, the first step in this derivation is to compute the retarded 
and advanced QD Green's functions. 

1. retarded and advanced QD Green's Functions 

One way to obtain the retarded QD Green's function (similar analysis can be done for 
the advanced QD Green's function) is by means of perturbation theory in the AC signal up 
to infinite order. For this purpose, one takes the QD Green's function ( |A12D as the bare 



QD Green's function. Using S^^(t) = —iVAccosuot (S^(^(t) = T.'^^it) — S^^(t), where 
S^^(t) = 0), as the retarded self-energy, the Dyson equation for the retarded QD Green's 
function reads 

/oo 
drgl^it - t) Vaccosuot gl^r - t') - 
-OO 

- dr dn gl^{t - t) Vaccosujot g'^^^ir - n) VaccosujoTi grf_^(ri - t') + ... (Bl) 

Z J —OO J —CO 



Eq. ([Bl| ) can be written as 

7T / VAccosciJoTrfT / VAccosuoTidri + ... . (B2) 

2 Jt Jt -' 

Performing the summation, one obtains the time dependent retarded QD Green's function: 



p=oo m=oo 

Y: E ■Jm{P)-UP)e-''^''e'"^''''slAt-t')- (B3) 

p=— oo m=— oo 

where Jm{P) is the Bessel function of order m and argument P = y^^. A similar procedure 
yields the time dependent advanced QD Green's function 

Gl,{t,t') = e-^^-"-o*e'^"""''*' gl^{t - t') = 

p=oo m=oo 

E E J™(/3)J,(/3)e"^P-«*e*'"-"*'g,%(t-t')- (B4) 

p=—oo m=—oo 



2. greater and lesser QD Green's functions 

The greater and lesser QD Green's functions are calculated by substituting the retarded 
and advanced QD Green's functions (Eqs. (|BB|), ( |B4| )) in Eq. (j^): 



./ 



(B5) 



p=— oo m=— oo L,R 

de gip<^ote-im<^oi'e-i<t-t')/^/^(e)r^/-f?(e) 

-oo 2n (e - pujo - eo,a + f EaeL,i?r„)(e - mcuo - eo,a - | EaeL.i? T, 

p=— oo m=— oo L,R 

J-oo 27r (e - pc^o - eO,a + f E^gL,/? Ta) (e - mUJo - eo,a - \ Ea&L,R Ta) ' 



p=oo m=oo 



p=—oo m=—oo L,R 

r^ de ^-ipu^ot^im^ot' ^-ie{t'-t) f^^^(^^^Y^/R 

7-00 27r (e - pcjo - eo,^ + f EagL,/? r„) (e - muo - eo,^ - | EaeL.i? T^) ' 
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p=—oom=—oo L,R 



fB8) 



D 27r (e - J9CJ0 - eo,a + f EaeL,i? rQ,))(e - mc^o - eo,a - f Eaei.ij ^a)) ' 
Once again we have taken eo,o- = eo,o- + Urid,^ being n^^^- the QD occupation. In the case of 
zero DC bias, the right and left Fermi-Dirac distribution functions are fifii^) = (e-Ep)/T 
where T is the temperature. The analytical expressions for the lesser and greater QD Green's 
function are obtained by integrating in the complex plane where the Fermi-Dirac distribution 
function can be written as a difference of two digamma functions which have poles in the 
lower and upper complex plane: 



1 1 



2 27ri 



IZ , , ,1 IZ 



k=oo 



/W - o = -7T- ^io + T^Z^) - ^(o - T^Z^) = E ^-T — - 1 , ,. , ^z ■ (B9) 



'2 2-nT' ^'2 2txT 



1 1 



h — ^^^ i -I- A- -I- 

fc=0 2 ' '^ 2-kT 2 ^ "^ ^ 2-kT 



The complex integrals are performed with the restriction t — t' > because we use them to 
evaluate a retarded quantity (the retarded self-energy in Eq. (12)). 
In order to abbreviate the notation, we define the following variables 

Op = -(eo,a + P^o - ^ Yl ^a)' ^™ = -(eo,a + muo + - ^ r„), 

^ oeL.R ^ a&L,R 

"^ " 2 "'' 2^""^' ^™ = 2 ^ 2^ "' ^ " "" 

Up = -(1 - tanh{—ap)), Vm = -(1 - tanh{—b^)). (BIO) 

With this notation, the analytic expressions for the lesser and greater QD Green's functions 



are: 



VAn Vac- , P—°° m=oo ApuJot ^-imoJot' 



G'>(t,t') = e-*^"""°V^^™"''* Y, J2 JU(3)Jpif3)- 



.Fi{ap,l,ap + l,C) Fi(7„,l,7^-hl,C) 



ftp 7" 



27ri (Mp - 1) e*"-(*-*')), (Bll) 



■Vac ■ . -Vac ■ ./ ^^^ '^=~ ipuot -imuot' 



p=-oom=-oo (P - "^)^0 - ^ EaeL.i? Tc 

Fi(a;,i,a; + i,c) Fi(7;;,i,7;; + i,c) 
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2711 Vm e-*"(*-*')), (B12) 



.VAn ■ -Va^ . , P~°° m=oo -ipLOot imujot' 



Fi(a;,i,a; + i,c) Fi(7;;,i,7;; + i,C) 



X 



-27rz(t;„-l)e-'''"(*-*')). (B14) 



"p /m 

The functions Fi are hypergeometric functions. Here 7* and a* are complex conjugates of 
7 and a respectively. Once the greater and the lesser QD Green's functions are obtained, 
the retarded QD interaction self-energy is given by: 

S2!f^(t,0 = -^UH{t - t')[G<At.t')Gl-At\t)G<At,t') - G>Jt,t')G<,{t',t)G>,{t,t')]. (B15) 

Here, we want to point out the non trivial dependence of the retarded self-energy on the 
parameters of the AC voltage. With this non trivial dependence the QD DOS strongly 
deviates from the usual single-particle Tien-Gordon behavior. [ p]| , P|] 
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FIG. 1. Complex time-contour. The times in the positive branch are r_|_ while times in the 
negative branch are r_. 
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t_ i\ L + 

FIG. 2. Self-energies of order U'^, S^^(t, t') and S^^(t, t'). The times in the causal branch 
are marked with a + symbol whereas the times in the anticausal branch are marked with a — 
symbol. Solid lines denote QD Green's functions in the Hartree approximation including coupling 
to the leads and AC potential. Dashed lines correspond to the on-site repulsion U electron-electron 
interaction. 
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FIG. 3. Time-averaged DOS in the strongly correlated regime, U = 2.57rr and T = O.OSr. (a) 
The AC frequency is ujq = r/4 ~ 2Tk, the solid line is the case in the absence of the AC potential, 
i.e., (3 = 0, where (3 = -^^. In this case the peak at Ep reaches a height of 0.93. The dotted line 
corresponds to /3 = 0.25 {Vac = 0.0625r ~ Tk/2). The dashed line shows the case of /3 = 0.5 
{Vac = 0.125r « Tk). The dot-dashed line corresponds to /3 = 0.75 {Vac = 0.1875r « 3Ti^/2). 
In the three cases the Kondo peak is slightly reduced by the AC signal. 

(b) ujq = r/2 ~ 4Tx. The solid line corresponds to (3 = 0, the dotted line to /? = 0.25 
{Vac = 0.125r « Tk). The dashed line shows the case of /3 = 0.5 {Vac = 0.25r « 2Tk) 
where the peak has been significantly reduced and the dot-dashed line shows the case of an intense 
signal where f3 = 0.75 {Vac = 0.375r ~ STk)- In this case the replicas of the Kondo become 
apparent (located at ujq and —too) and the peak at Ep has been strongly reduced. Figure (c) 
corresponds to ujq = T: the solid line is for the case (3 = 0; the dotted line shows the case /3 = 0.25 
{Vac = 0.25r ~ 2Tk)- At this AC intensity the first satellites of the Kondo resonance show up. 
The dashed line corresponds to (3 = 0.5 {Vac = 0.5r ~ 4Tx), in this case the time-averaged DOS 
at Ep has been suppressed below 0.5. The dot-dashed line corresponds to a very intense signal, 
P = 0.75 {Vac = 0.75T « 6Tk) , the Kondo peak in this case has completely disappeared. Insets: 
Blow up of the time-averaged DOS near w = 0. 
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FIG. 4. (a) Conductance as a function of Vac ™- the correlated regime where, U = 2.57rr and 
T = O.OSr for four AC frequencies. 
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FIG. 5. (a) Conductance as a function of ujq in the strong correlated regime U = 2.57rr and 
at T = O.OSr. The solid line shows the case of Vac = 0.25r ~ 2Tk, the dotted line corresponds 
to the case of stronger AC intensity Vac = O.SP ~ 4Tk (b) Time-averaged DOS for the cases 
of ujQ = 0.375r (solid line), loq = T/2 (dotted line), ujo = 0.625r (dashed line) for a fixed AC 
amplitude Vac = O.SP. 
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FIG. 6. (a) Conductance as a function of cjq in the Kondo regime, U = 2.57rr and T = 0.05r. 
The sohd hne shows the case of Vac = T, the dotted hue corresponds to the case of Vac = l-SF 
and finally dashed line shows the strongest AC intensity Vac = 2r. 

(b) Time- averaged DOS for the strongest field amplitude Vac = 2r for different large AC frequen- 
cies. The solid line shows the case of ujq = 3.5r (/? = 0.57) where the Kondo peak has the lower 
height, the dotted line shows the case loq = AT {(3 = 0.5), the dashed line corresponds to loq = f>T 
{(3 = 0.22) and the long-dashed line depicts the case of ljq = GT (/? = 0.16). 
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FIG. 7. Time-averaged DOS at Ep for U = 2.5-7rr as a function of temperature. The solid 
line shows the case in the absence of the AC potential. The dotted line corresponds to a case 
of a low AC intensity Vac = 0.25r ~ 2Tk (/3 = 0.25) and the dashed line shows a case with 
higher AC intensity Vac = O.SF ~ 4Tk (/3 = 0.5). Both of them correspond to an AC frequency 
uJo = T » Tk- 
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